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Abstract 

We demonstrate an efficient nonlinear process in which Stokes and anti-Stokes 
components are generated spontaneously in a Raman-like, near resonant me- 
dia driven by low power counter-propagating fields. Oscillation of this kind 
does not require optical cavity and can be viewed as a spontaneous formation 

of atomic coherence grating. 
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Theoretical and experimental work of past few years on atomic coherence and interfer- 
ence has demonstrated a potential to improve significantly the existing nonlinear optical 
techniques [jl|. In the present Letter we demonstrate efficient parametric generation ac- 
companying the spontaneous formation of the coherent superposition states. Specifically, 
the present work reports the observation of the spontaneous parametric self-oscillation in 
resonant Raman media. Such generation does not involve optical cavity and appears under 
remarkably simple circumstances, when two low-power counter-propagating fields interact 
with the medium. Oscillation manifests itself as Stokes and anti-Stokes components gener- 
ated with a frequency shift corresponding to that of the Raman transition. As the oscillator 
goes over threshold, dramatic increase and narrowing of the beat note between the input 
field and generated components takes place. 

The principal possibility of mirrorless parametric oscillation with counter-propagating 
signal and idler fields has been suggested in 1960's by Harris 0. The original proposal 
based on non-degenerate frequency mixing has not been realized up to now due to small 
values of nonlinearities in available materials and difficulties in achieving phase matching . 
It is easier to achieve mirrorless oscillation in degenerate four-wave mixing. The possibility 
of self-oscillation in such interactions has been predicted in [^], and a number of the related 
effects, such as conical emission or transverse pattern formation have been observed in a 
vapor driven by very strong, off-resonant counter-propagating laser beams |^. Workers in 
the field have also noted the importance of Raman nonlinearities in the early experiments 
on polarizations instabilities 0. 

As compared to the above work the presently reported results utilize atomic coher- 
ence gratings [^|[ in a resonant double-A atomic system (Fig. la). Coupling of counter- 
propagating Stokes and anti-Stokes fields via such a grating appears to be the main physical 
mechanism resulting in Raman self-oscillation p. Similar to several related studies p|,|7|,p|JT0[] 
the present work operates in a so-called strong coupling regime in which nonlinearities can 
not be derived from a usual perturbation expansion. In this regime quantum coherence 
and interference have a profound infiuence on nonlinear parametric amplification. For ex- 
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ample, linear and non-linear absorption of parametrically generated fields can be controlled 
and the phase mismatch, inherent in all non-degenerate parametric interactions involving 
counter-propagating fields, can be easily compensated by a large dispersion accompanying 
resonances in phase-coherent media. It is important that due to quantum interference the 
strong-coupling regime was reached in the present work with a very low driving power. The 
present results are therefore directly related to recent theoretical studies on few photon quan- 
tum control fill], switching |jl2|, and quantum noise correlation |jl3|, and can potentially be 



used to study interactions of a very low-energy fields and for suppression of quantum noise. 
Moreover, if extended to resonant molecular vapor the present approach might be useful 
for efficient Raman frequency shifting. Likewise, narrow-linewidth signals may also be of 
interest e.g. for optical magnetometry. 

The present experiments follows our previous work on atomic coherence effects in opti- 
cally driven A systems in Rb 0,0. Studying the detailed lineshape of these signals, we 



found that under certain conditions the Raman amplification can actually turn into a 
coherent self-oscillation of the Stokes and anti-Stokes components. An essential element for 
the oscillation to appear is the existence of the two driving fields {Ef,Ei,) propagating in 
the opposite directions. We first discovered this oscillation as the result of simple Fresnel 
refiection from the rear window of the Rb cell. Figure lb shows the simplest experimen- 
tal configuration that produced Raman oscillation. A beam from an extended-cavity diode 
laser passes successively through an optical isolator (I), a focusing lens, a heated Rb cell, a 
partially refiecting mirror and onto a fast photodiode (PD). The signal from the photodiode 
is detected using a microwave spectrum analyzer (SA). The partially refiecting mirror (M) 
is used to retro-refiect some of the transmitted beam back through the Rb vapor. With 
proper tuning of laser frequency, the backward beam causes the Raman gain peak to grow 
to the point of oscillation threshold. When the self-oscillation occurs the detected Raman 
beatnote signal at a frequency of hyperfine splitting {ujhfs) increases in amplitude by as 
much as 60dB and its linewidth narrows from 200 kHz to less than 300 Hz (Fig. 2a). Under 
appropriate conditions the beatnote linewidths as narrow as 100 Hz FWHM were observed 



(Fig. 3a). This is much narrower than the usual broadening mechanisms for Raman transi- 
tions under the present conditions conditions (primarily, transit broadening ~ 50 kHz 
and power broadening ~ 500 kHz). The oscillation occurs without any cavity enclosing the 
cell. We have been careful to eliminate possible extraneous sources of feedback to lasers or 
other optical and electronic elements. 

In order to study the physical origin of the oscillation process we carried out a series of 
experiments, where instead of reflecting the incident laser light we injected laser beams with 
different frequencies from the opposite directions (Fig.lc). We found that the oscillation 
occurs readily if the forward and backward fields are tuned to the different ground state 
hyperfine level as diagramed in Fig. la. It is more difficult to make the system oscillate 
if the backward beam is tuned to the same frequency as the forward beam. If tuned to 
different ground state hyperfine levels, the oscillation was observed with the backward beam 
coupling either the same (-P1/2) or different (P3/2) upper-state fine-structure levels as the 
forward beam. In our two-laser experiments it is was easy to see oscillation for both ^^Rb 
and '^''Rb isotopes. When oscillating, the Rb vapor can convert as much as 4% of the total 
input power into the frequency shifted Stokes and anti-Stokes sidebands. 

Typical conditions to observe the oscillation with a single laser beam are: ECDL tuned 
in the wing of the Doppler profile of the ^^Rb, Dl line (F=2 to F'=3) transition, ~ 10 mW 
of power in spot size about ~ 500 microns, 5 cm long Rb cell operated at 75-95 C, and 
between 10 and 80 % of the driving power is retro-refiected back through the Rb cell. In 
the two-laser experiments oscillation was observed for driving input powers 2 — 10 mW, spot 
sizes 0.1 — 2mm, and cell temperatures 65 — 100°C. 

The oscillation frequency shift {ujq) does change somewhat with laser tuning (typical case 
was 30 Hz per MHz of laser tuning) and with angle between the forward and backward beams 
p!5| . However, the oscillation frequency always remains within the bandwidth (few hundred 
kHz) of the power broadened and shifted single-beam Raman gain peak. The oscillation 
prefers, but does not necessarily require a circularly polarized beam, and the gain is largest 
with zero applied magnetic field. 
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We have analyzed the characteristics of the forward and backward beams by making 
beatnotes with independently tuned laser sources, and by using optical cavities to analyze 
the spectra. We found that the field components at frequencies of the forward and backward 
driving fields {uf or z/fe) are surrounded by generated first order Stokes and anti-Stokes fields 
at frequencies Uf h ± uq. In certain cases second order components have been seen as well. 
The generated components produce, in general, asymmetric spectrum. In particular, in 
cases when forward driving field is tuned to e.g. upper ground state hyperfine sublevel and 
backward driving beam is tuned to the lower hyperfine sublevel, the anti-Stokes component 
observed in a forward direction is much more (~ 20 — 30dB) intense than the Stokes one. 

These observations suggest that the actual oscillation mechanism is somewhat different 
from (although related to) that studied theoretically in Ref . . That work involved only one 
pair of counter-propagating components. Motivated by experimental results, we consider a 
theoretical model in which atoms in a double A-type configuration are interacting with six 
optical fields. These include two counter-propagating driving fields with frequencies up, 
and complex slowly varying amplitudes Sp and Sb', anti-Stokes and Stokes components with 
frequencies i^i^s = up ±ujq propagating in the forward direction {81,83), and corresponding 
components with frequencies 1/2,4 = '^b ± t^o propagating in the backward direction {82,84). 
The field is then written as E = Y,i{£i^~^^'^^^~^''^^^ + c.c.)/2. Below we focus on the linear 
theory describing the oscillation threshold. Hence, all generated components are treated 
to first order only and saturation effects are disregarded. These assumptions allow us to 
truncate the infinite hierarchy of equations. The resulting polarization can be written in 
the form P = Ei(^ie-^('^'*+^»^) + c.c.)/2. We are interested here in polarizations at the 
Stokes and anti-Stokes frequencies, which are related to the field components by 4 x 4 sus- 
ceptibility matrix Xmn- Vm = eoXmnew{ikrnnr)8n, where P = [Vi, V2, (P3)*, {'P^Yf, 
8 = [81, 82, {83)*, {84)*]^, and kij are representing free-space wave vector mismatch. 
For the present problem the matrix elements of [x] were calculated explicitly for each veloc- 
ity group and averaged over Maxwellian velocity distribution. In the present calculations we 
consider fields interacting in a slab of medium of the length L. Assuming that the solution 
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is homogeneous in transverse directions leads to {kij)± = 0, and the evolution along the 
longitudinal direction z is described by: -^Si = i/{2eo){ki)zVi. The appropriate boundary 
conditions are taken to include a weak "seed" input (S) at anti-Stokes frequencies (corre- 
sponding to e.g. spontaneous emission, or vacuum field). 

Before proceeding with comparison of experiment and theory we illustrate the origin of 
the oscillation. To this end, let us assume that absorption of the driving fields is negligible, 
and there is no inhomogeneous broadening. Furthermore, we disregard the coupling of the 
forward (backward) driving field with all transitions except for c ^ a (6 — a') and assume 
that the detuning of the backward driving field from respective single photon resonance (A^) 
is much larger that the corresponding detuning of the forward drive. In such a situation 
only forward anti-Stokes (Si) and backward Stokes (84) fields are involved into nonlinear 
interaction (Fig.4a, 0). In this case: d£i/dz = a^Bj^ with = {1,4} and an = 

-Villbc + KuJo - oJhfs) + i{\^B? - \VLf\^)/ ^b\/\^f\'^ - iku, au = ir]i[QB^F/{^B\^F\'^)], 

= irii[fl*^fl*p/{AB\^F\'^)], (^u = —ik^i. Here rji = (A;j)23/(87r^)A^(Aj)'^7i, where Aj is a 
wavelength of the ith field component and 7^ is the radiative decay rate on the transition 
coupled by this component of the field. is atomic density, and Qf,b are Rabi- frequencies. 

When the phase matching condition is satisfied (Im((5a) = 0, 6a = (an — au)/2), we 
find: 

S,{L) ~ S,{0r ~ / — -, (1) 

da sm(sL) — s cos(sL) 



where s = \J aua4^i — (5a)^, and the unimportant proportionality constants have absolute 
values of the order of unity. These solutions diverge if tan(sL) = s/(5a), which indicates the 
onset of mirrorless oscillations. Note that the latter condition can be satisfied if rii\VLpVLB\ > 
rji'ybcl^Bl, which is identical to a strong coupling condition of Refs. |[l|J^-[T^. Let us examine 
now the phase matching. Close to the two-photon resonance we have: 

k{ujo - uJhfs - + + h-ki- ki]^ = 0, (2) 

where ^ = (|f2j7p— \Qb\'^)/^b is a function of drive power and detuings representing a 
light shift, and n = crii{ki)z/\flF\'^ ■ It is interesting that this equation resembles closely the 
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frequency pulling equation of the usual laser theory, with frequency stabilization coefficient 
K. The first term in the left-hand side corresponds to atomic dispersion, and the second 
describes the geometrical phase mismatch. This contribution is proportional to the Raman 
transition frequency (see inset to Fig.4) and also depends on relative angles between driving 
beams. Hence it plays a role analogous to the cavity shift. Note, however, that under the 
typical oscillation conditions stabilization coefficient k ~ c/(7fcc-^) ^ 1 and the oscillation 
frequency is locked to the light-shifted Raman transitions frequency. 

This implies that in the example considered above the physics behind the oscillation 
phenomenon is the coupling of the counter-propagating Stokes and anti-Stokes fields via 
spontaneously created atomic coherence p^c- On one hand such coupling results in scattering 
of these fields into each other, thus forming an effective feedback. One the other hand, this 
process is also accompanied by parametric amplification. When both of the effects are 
present self-sustained oscillation can occur. 

Let us consider the influence of other nonlinear processes on mirrorless oscillation. When 
only a forward driving fleld is present nonlinear interaction results in the coupling between 
forward (or backward) propagating Stokes and anti-Stokes fields (Fig. 4b) leading to coher- 
ent Raman scattering and amplification of the co-propagating pair of fields in the vicinity 
of two-photon resonance Oscillation is not possible in this case, since no effective 

feedback is present. However, when coherent Raman scattering exists in addition to the 
coupling between counter-propagating Stokes and anti-Stokes components, it can result in 
lowering the oscillation threshold. The process shown in Fig. 4c represents a different type of 
parametric interaction. It leads to the scattering of the counter-propagating anti-Stokes (or 
Stokes) waves into each other, which does not change the total photon number of weak fields. 
Consequently, it alone can never lead to the oscillations. However, oscillations can emerge 
if in addition to the parametric energy exchange additional amplification mechanisms (e.g. 
coherent Raman scattering) are present. 

In general, for the detailed comparison of the theory and experiment all of the six pro- 
cesses of the type shown in Fig.4 should be taken into account. They give rise to simultaneous 
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generation of all components in both directions. To make a comparison we have solved the 
full system of propagation equations numerically, taking into account Doppler broadening, 
and propagation of all fields. The results (Fig. 2b) show good qualitative agreement with 
experiments. 
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L. HoUberg and V. Velichansky. The authors warmly thank them as well as M. Fleischhauer, 
P. Hemmer, S. Harris, A. Matsko, V. Sautenkov, and G. Welch for useful discussions, and 
T. Zibrova for valuable assistance. We greatfuUy acknowledge the support from the Office 
of Naval Research, the National Science Foundation, and the Air Force Office of Scientific 
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FIG. 1. (a) A prototype 4- level model for self-oscillations. In general, we assume that each 
driving fields couples both of the ground states. The upper levels of this double-A system can 
represent some manifolds of states. In a particular case, a and a' can also represent an identical 
state. In the experiment states c and h are hyperfine sublevels of the Rb ground state 5>S'i/2- (b,c) 
Experimental setups (schematic). 
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FIG. 2. a) A typical signal recorded by a fast photo-diode. Curve a is recorded with only 
forward driving beam present. Self-oscillations occur in the presence of the forward and backward 
driving fields (curve b). Parameters are: cell temperature 92 ° C; forward driving beam with power 
10 mW and spot size 1.5 mm is detuned by 800 MHz to the red side of F = 3 ^ F' = 3' transition 
of Di line; backward driving beam with power 2.5 mW and spot size 1.5 mm is detuned by 2 GHz 
to the blue side of F = 2 ^ F" = 3" transition of D2 line, b) Calculated signals corresponding to 
the experimental conditions of Fig. 2a. (At the point corresponding to parametric oscillation linear 
theory predicts infinite growth). 
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FIG. 3. a) A typical beat signal at 3.034 GHz recorded when the frequency of the driving laser 
is locked to a reference cavity, b) A typical dependence of the generated anti-Stokes power as a 
function of the input driving power in a vicinity of oscillation threshold. 




FIG. 4. Examples of different types of nonlinear processes contributing to self-oscillations, a) 
Direct nonlinear coupling between counter-propagating Stokes and anti-Stokes fields; b) Coherent 
Raman scattering; c) Parametric energy exchange between counter-propagating anti-Stokes fields. 
For each of the processes (a-c) there exists a complimentary process of the same type involving 
other pair of weak fields (e.g. £^2,3 in Fig.4a). Inset illustrates wave vector mismatch for the process 
(a). If all fields are propagating along the z axes [/cq^ + fc^ — fci — fc4]2 = 2u}q/c. 
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